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Abstract 

In film growth on quasicrystalline surfaces, the epitaxy-imposed ordering cannot 
compete with the stable bulk phases of thick films due to absence of translational 
order in the structure of the substrate. Energetically, this renders the formation 
of crystalline domains in the native structure of the film material more favorable 
while their global orientation is prescribed by the quasicrystalline order. We 
present experimental results on the dissociative chemisorption of oxygen on the 
decagonal surface of AlyoCoisNiis as well as molecular-dynamics simulations of 
the diffusion of adatoms on the surface of the partially covered substrate. 
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1. Introduction 

In hcteroepitaxial growth, the properties of the film strongly depend on 
the amount of structural mismatch with the substrate as well as on the film 
thickness. The latter adversely infiucnccs the elastic energy of the strained film 
giving rise to a three-dimensional growth regime and the formation of islands 
with limited sizes. This phenomenon can be utilized to fabricate nanometer- 
size domains by depositing material on a substrate with a large lattice mismatch 
of typically in excess of 5 % An extreme case is encountered by growing 
crystalline material on a quasicrystalline substrate for which no translational 
periodicity exists [2j . We will outline below our experimental and computational 
efforts to elucidate the processes encountered in such an extreme case of lattice- 
mismatch in heteroepitaxy. 

The decagonal quasicrystal AlypCoisNiis (d-AlCoNi) has a structure which 
is quasiperiodic in two dimensions, periodic in the third dimension, and dis- 
plays decagonal diffraction symmetry. Investigations of its bulk structure have 
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revealed a columnar prismatic morphology with the column axis parallel to the 
periodic direction, the tenfold-symmetry axis OH]. The decagonal surface is 
perpendicular to the tenfold-symmetry axis and is very suitable as a substrate 
for epitaxial growth due to its flatness |5Jin]. 

We have previously reported that Al adsorption on the tenfold-symmetry 
surface results in 3-nm large Al(lll) domains for coverages in excess of a few 
monolayers (ML). The azimuthal orientation of the domains is governed exclu- 
sively by the aperiodic symmetry of the substrate [7] . We also have studied nu- 
merically the effect of the relative strength of mutual interactions of adsorbed 
atoms (adatoms) with respect to their interactions with the substrate atoms 
[SJ [H] on the structure of the growing adsorbate layers and how the structure of 
the adsorbate layer (adlayer) depends on the rate at which adatoms thermally 
equilibrate with the substrate |10j . Deposition is a non-equilibrium process and 
the ratio of the thermal relaxation to the deposition rate of the adatoms has 
influence on the resulting morphology of the adsorbate [TT] . We found in par- 
ticular that the thermal relaxation rate and relative strength of the interaction 
of adsorbate atoms with each other and the substrate was able to change the 
surface growth mode from cluster to layer-by-layer type. We attributed this to 
(i) the extent to which arriving adatoms can diffuse on the surface before being 
thermalized and (ii) how strongly already present adatoms can steer this dif- 
fusion process. Our purpose in this report is to further elucidate this behavior 
by focusing on the diffusion of deposited adatoms on the substrate-adsorbate 
surface. We also present experimental results on the initial stages of the high- 
temperature oxidation of the tenfold-symmetry surface of d-AlCoNi resulting 
in the formation of hexagonal domains of dissociated oxygen. The thickness of 
the oxygen adlayer as well as its properties are determined by the conditions at 
which oxygen interacts with the surface. Hence, this is a self-regulating system 
for which the physics and chemistry at the interface determine the flnal surface 
constellation. 

2. Experimental 

In this report we have used a d-AlCoNi quasicrystal with dimensions 5 x 3 x 
1 mm"^ which was mounted on a resistive heater, while the sample temperature 
was measured with a chromel-alumel thermocouple (K-type) pressed onto the 
sample surface |12j . The quasicrystal surface was cleaned in ultrahigh vacuum 
by cycles of sputtering with Ar+ ions (1.5 keV, 5x10^^ A/mm^) at 670 K 
and heat treatment at 900 K for 30 minutes. The sample preparation was 
monitored by the quality of the low-energy electron diffraction (LEED) pattern 
and by inspecting the scans of Auger electron spectroscopy (AES). A three-grid, 
back- view, display- type LEED system, operating at low microampere primary 
currents, had a total opening angle of 102±2° |13j . Thus, a momentum transfer 
of 2.83 ± 0.04 could be detected for 50 eV electrons at normal incidence. 

The diffraction patterns were recorded by a 16-bit charge-coupled device 
camera after the sample was cooled down to room temperature. The position 
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Figure 1: LEED patterns at a primary-electron energy of 50 eV obtained from (a) the clean 
and (b) the exposed tenfold-symmetry surface of d-AlCoNi to 100 Langmuirs of oxygen at 870 
K. The patterns are shown after normalization by the response function of the display and 
recording system in order to eliminate spurious contributions to the image. 

and the size of the Bragg spots in LEED observations are used to extract real- 
space information about the atomic structure, the size, and the orientation 
of the surface textures. The AES spectra were recorded using a cyHndrical 
mirror analyzer operating with a constant relative resolution of /S.E/E — 0.8%. 
Experimental details have previously been reported [HI [71 [T51 [TU] . 

3. Oxygen adsorption at high temperatures 

Fig. shows a LEED pattern from the clean tenfold-symmetry surface of 
d-AlCoNi at 50 eV. The pattern displays tenfold symmetry with spot profiles 
having a total width of about 0.5° corresponding to an average terrace size of at 
least 150 A. This value is comparable with the coherence length of the primary 
electrons used here and, therefore, is the upper limit of the domain size that 
can be determined in our apparatus |13j . 

The LEED pattern depicted in Fig. [TJd is obtained from the surface af- 
ter exposure to oxygen partial pressure of 1 x 10~* mbar at 870 K for 1000 
s. We discern thirty new diffraction features placed on a polar circle of about 
40.5 ± 0.5°. In accordance with the local symmetry of the quasicrystalhne sur- 
face, these spots represent a sixfold-symmetric pattern, repeated five times in 
equal azimuthal increments of 27r/5. Each sixfold-symmetric pattern is due to 
a hexagonal mesh with an interatomic distance of 3.08 ± 0.03 A, as extracted 
from LEED. We note that his value is 10% larger than the interatomic 0-0 dis- 
tance in different phases of AI2O3 [TB]. A lattice expansion of a similar amount 
has also been encountered in oxide films grown on Al7oPd2oMnio |T5]. It is 
important to note that the intensity of the new diffraction features do not show 
any polar smearing compared to those from the quasicrystal, while they are 
extended azimuthally by 2 — 3° . We may interpret the polar precision as an un- 
expectedly large domain size of about 150 A, while the azimuthal spread arises 
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Figure 2: Auger spectra obtained from the clean (red, bottom) and oxygen-adsorbed (blue, 
top) tenfold-symmetry surface of d-AlCoNi in the energy region of the Al as well as O, Co, 
and Ni transitions. Spectra for the high-energy region are enlarged by a factor of 3. 

as a consequence of the relaxation of the crystaUine adlayer compared to the 
quasicrystaUine substrate. The latter facilitates the growth of the adlayer in this 
extraordinarily large size. We note that the diffraction spots characteristic of 
the clean quasicrystaUine surface are still observable indicating that the surface 
film is rather thin (few A), even thinner than 5 A, the thickness of oxide films 
grown on the fivefold-symmetry surface of Al7oPd2oMnio [H]. The observation 
of these spots at the same place further indicates that the formation of the sur- 
face layer preserves the quasicrystaUine order at the interface. The azimuthal 
alignment of the spots with respect to the substrate is shifted by 6° compared 
to that observed for Al films deposited onto the same decagonal surface [7] . 

Fig. [2] shows Auger spectra near the Al 2pVV, O IsVV, Co and Ni 2p3d3d 
transitions for the clean (bottom) and oxygen-exposed (top) surface at 2.4 keV 
excitation energy. The spectrum obtained from the clean surface is dominated 
by the Al transition at around 70 eV similar to that obtained from clean Al. 
For an oxidized Al sample this peak is found shifted to lower energies by almost 
10 eV due to a strong electron transfer from Al to O [17]. In the present case, 
however, no energy shift for the Al signal is observed after oxygen adsorption, 
but some reduction of the intensity. This observation indicates that there is no 
appreciable electron transfer from the metal site to oxygen, i.e., oxygen is in 
the chemisorbed state. The tiny shoulder at around 60 eV in the top spectrum 
may signal the presence of a minute amount of aluminum oxide and hence the 
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early stage of oxidation. Around 500 eV the spectrum from the oxygen-exposed 
surface shows the oxygen IsF^ signal. Both surfaces display 2p3d3d transitions 
for Co and Ni at higher kinetic energies. Except for a slight intensity reduction, 
no chemical effects can be discerned in this energy region for the transition-metal 
signals indicating that oxygen binding to Co or Ni is negligible. 

Taking into account the features of the LEED pattern and the results of 
AES, the adsorbed oxygen may be just one single ML thick. It is interesting 
to note that while a ML of noble gas forms a pseudomorphic layer [18] , oxygen 
adsorption distinctly forms a hexagonal lattice on the quasicrystalline surface. 

4. Simulation of adatom diffusion on the adsorbate-substrate surface 

We perform molecular-dynamics simulations of how deposited adatoms dif- 
fuse on the adsorbate-substrate surface in order to understand the morphology 
of the growing adlayer. We use a 6 000-atoms substrate model of the decagonal 
AlyoCoisNiis quasicrystal with dimensions 150 x 150 x 2.05 [101 EH]. The 
large lateral size of the unit cell reduces artefacts generated by boundary condi- 
tions. For the dimensions given above, 1 ML of adatoms corresponds to about 
3100 particles. In the simulations we inject adatoms in batches of 180 particles 
onto the substrate. The quasicrystalline bilayer is kept rigid, but adatoms axe 
allowed to interact with each other and with the substrate atoms. The relative 
strength of the adatom-adatom and adatom-substrate interactions is controlled 
by a dimcnsionless parameter rj. Experimentally, it is found that the adsorp- 
tion of Al on AlCoNi yields qualitatively similar results as the adsorption of 
nobel gases such as Xe on the same type of substrates [7l [18] . These findings 
imply: (i) given that Al and Xe have very different chemical binding properties, 
the mechanism of domain formation does not depend strongly on the details 
of the interaction, which motivates our use of a simple Lennard- Jones interac- 
tion potential (ii), the interaction between adsorbate Al atoms and Al in the 
substrate might be different from the Al interactions in crystalline Al and thus 
the strengths of the interaction of Al adatoms with atoms in the adsorbate and 
substrate are not necessarily equal. 

In real experiments, hot impinging adatoms loose their excess kinetic energy 
to the substrate and eventually reach thermal equilibrium. In our simulations, 
we account for this process by introducing damping and treating the damping 
rate 7 as a parameter by which we can change the thermal relaxation times. 
Numbers given below are in units in which lengths are measured in angstrom 
and the mass of the adatoms is set to one. 

Numerical details of the implementation are as follows. We let Adatom 
particles i interact with each other and the atoms a of the rigid quasicrystalline 
substrate via a Lennard Jones potential. 



with parameter values of cr = 2.55 (in accordance with the nearest neighbor 
distance of Al atoms in the fee structure) and e = 0.25, which establishes an 
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energy scale (eV). The equations of motion of the adatoms are given by. 



^ VF(|ri ~ r„|) - r;^ VV^drj - ri|) + f,(T), 



(2) 




where the parameter 77 is the relative strength of the adatom interactions with 
respect to the substrate and (T) is the damping term for the thermal relaxation 



with 7 being the thermal damping rate, T the temperature measured in units of 
e and f being the velocity unit vector of the particle [H] . We do not distinguish 
between the different atomic types in the substrate. The equations of motions 
are integrated using an adaptive step size fourth order Runge-Kutta scheme, 
[22]. Given the relaxational nature of the particle dynamics, the choice of an 
adaptive steps size controller turned out to be an efficient way of maintaining 
numerical precision at a minimum number of integration steps. 

We have already shown that by varying 7 and 77, different surface growth 
morphologies including layer-by-layer and cluster-type growth can be observed 
[3]. We attributed this finding to the fact that a higher damping rate, which 
corresponds to a short thermal relaxation time, causes the particles to hit the 
substrate and stick, very similar to rain drops wetting a surface, while the 
effect of varying ij determines how strongly already deposited adatoms can steer 
impinging adatoms during their diffusive motion |10j . We were able to explain 
how deposition at high 7 along with low 77 values, e.g., 7 — 0.25 and 77 = 1, gives 
rise to layer-by-layer type of growth, while keeping 7 constant and increasing 
77 will increase the probability of arriving adatoms to be deposited on already 
present ones giving rise to a cluster-type growth as observed in our simulations. 
We also found that well-ordered domains can grow at lower r/ values if the 
damping rate is reduced, so that impinging adatoms can diffuse longer on the 
substrate surface and, moreover, with increasing coverage, there will be more 
accessible favorable binding sites, since they will include both already deposited 
adatoms as well as substrate sites that have remained exposed [TU]. Fig. |3] 
demonstrate the formation of ordered and orientationally locked local hexagonal 
domains. We have performed a Voronoi tesselation of the the third adlayer above 
the interface and color-coded the Voronoi hexagons associated with sixfold- 
coordinated atoms that show close to perfect local order also in the bond-angles. 
The color of each hexagonal Voronoi cell is determined by the angle relative 
to the X-axis and corresponds to a range from to 60°. Multiple domains 
with relative orientation of 6°, 18°, 30°, and 42° are readily discerned. The 
fixed orientation angles of the domains are a result of the registry with the 
underlying tenfold symmetry of the quasicrystalhne substrate. An excerpt of 
the three-dimensional configuration of the adsorbate film is shown in Fig. [5] 

Here we aim to further substantiate these observations. Specifically, we 
analyze the deposition and the subsequent lateral diffusion of an additional set 
of 180 adatoms (0.06 ML) injected onto the adsorbate surfaces when 0.5, 1.5, 
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Figure 3: Voronoi tessalation of the third adlayer, as counted from the quasicrystal interface 
at 3.5 ML coverage, low damping rate 7 = 0.0625, and strong adatom interactions rj = 2. 
Voronoi cells associated with atoms exhibiting close to perfect local sixfold coordination have 
been color-coded according to the angle that the corresponding hexagons make with the x-axis. 
Notice the formation of nano-scale well-ordered domains with fixed orientations of 6°, 18°, 30°, 
and 42° . The fixed orientational angles are a direct consequence of being in registry with the 
tenfold symmetry of quasicrystalline substrate. Another possible domain orientation of 56° is 
absent. 

and 2.5 ML of adatoms have already been deposited. We ask how the lateral 
diffusion of impinging adatoms depends on the already present coverage, the 
damping rate 7, and the relative strength of the adatom interactions rj. We track 
the lateral diffusion of the impinging adatoms over a time scale sufficiently long 
to reach thermal equilibrium with the substrate and measure the displacement 
from the initial lateral positions as a function of time. 

Fig.|4]shows the root mean squared lateral displacement Ar^^^s averaged over 
impinging adatom trajectories that stick to the substrate surface as a function 
of time for rj — 2 and 7 values of 0.0625, 0.125, and 0.25 depicted as solid, dotted 
and dashed lines, respectively. The dependence of the diffusion on the coverage 
of already deposited atoms on the surface is shown, color-coded as follows: 0.5 
ML (red), 1.5 ML (blue), and 2.5 ML (green). We clearly see that the maximal 
diffusion distance decreases sharply with increasing thermal relaxation rate, as 
one expects. On the other hand, for a low relaxation rate (7 = 0.0625, solid 
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Figure 4: Root-mean-squared lateral diffusion distance against time. Shown are curves from 
three different damping rates 7 and and surface coverages: 7 = 0.0625 (solid curves), 7 = 0.125 
(dashed curves), and 7 = 0.25 (dotted curves), for coverages of 0.5 (red), 1.5 (blue), and 2.5 
(green) ML. See text for further details. 



lines), the diffusion distance drops sharply from 40 A at a coverage of 0.5 ML to 
8 A at 1.5 ML along with a lesser drop to 6 A at 2.5 ML, while the dependence 
of the diffusion distance on the substrate coverage becomes less pronounced for 
higher damping rates. 

It turns out that for 77 — 1.2 almost the same diffusion is obtained at 1.5 
and 2.5 ML and coverages, while the behavior at 0.5 ML coverage yields shorter 
distances. Thus, in general, it seems that the effect of 77 on the diffusion range 
is less pronounced or present only at low coverages. Nevertheless, we know from 
our previous work that the adsorbate covered surfaces reveal different morpholo- 
gies. In Fig. [5] we show 150 x 50 x 15 A'^ sections of the adsorbate for values 
of 7 = 0.0625 (left column) 7 = 0.25 (right column) and 77 = 2 (top row) and 
77 = 1.2 (bottom row) after deposition of 2.5 ML. Although we find compara- 
ble diffusion distances for configurations with different 77 (configurations in the 
same column), the morphology of the corresponding surfaces is different. For 
both 7 — 0.25 and 0.0625, the surfaces with stronger relative adatom interac- 
tions, rj — 2 exhibit enhanced well-ordered hexagonal domains, which are less 
pronounced for 7; = 1.2 (bottom row). The adlayer formed under low damping 
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Figure 5: Structure of the adsorbate surface at 2.5 ML coverage at a damping rate of 7 = 0.25 
at different strength of the adatoms interactions relative to the substrate. The top bottom 
panels show the resulting adlayer morphology for rj = 2 and 1.2, respectively. Notice the more 
locally ordered structure for stronger adatom interactions (ry = 2). Also, for rj = 2 the adlayer 
is more corrugated. 

and strong adatom interaction shows the largest amount of corrugation along 
with local hexagonal order (top left) . 

5. Discussion and Conclusions 

The dissociative adsorption of oxygen has already been observed on Al(lll) 
|20j , where the formation of a hexagonal oxygen structure is compatible with the 
substrate symmetry. In the present case, however, the substrate, while contain- 
ing Al, has a quasicrystalline structure on which the sixfold symmetric oxygen 
adlayer assembles. While the Al-0 interface for Al(lll) is stable (nonreactive) 
up to temperatures of 470 K, Al in the quasicrystalline Al-Co-Ni matrix remains 
nonreactive at temperatures close to 900 K. 

We observe the formation of unusually large oxygen islands on the qua- 
sicrystalline substrate. The influence of the latter, beyond providing the correct 
physical and chemical environment for the formation of the novel oxygen is- 
lands, consists of aligning these islands in a quasicrystalline matrix. Each island 
is quite large and low-energy electrons used in imaging the surface do not show 
a coherent interaction between the islands. 

In the simulations, we have used a quasicrystalline bilayer perpendicular to 
the periodic direction of the structural model to investigate the lateral diffusion 
of adatoms as they search for energetically favorable locations. We find that the 
morphology of the developing adlayer can be explained in terms of the thermal 
relaxation times for the impinging adatoms to reach thermal equilibrium and the 
relative strength of the adatom interactions with each other and the substrate. 
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Short thermal relaxation times generally lead to a layer by layer growth mecha- 
nism, in which the randomly incident adatoms hit the substrate and diffuse very 
little. The diffusion range is also found to diminish with the coverage, but this 
behavior is less pronounced at small thermal relaxation times. For coverages 
of 1.5 ML and more, the relative strength of the adatom interaction seems to 
have little effect on the diffusion range. An inspection of the forming adlayer 
however reveals that its morphology depends on how strongly adatoms interact 
with each other relative to the quasicrystaline substrate atoms. As the relative 
interaction strength of the adatoms is increased three dimensional cluster-type 
growth and island formation is observed. As we showed before, the structure 
of the adlayer either is an imperfect continuation of the quasicrystalline sub- 
strate maintaining local five fold coordination (QC-like phase), or it consists of 
nano-scale fee crystallites with their (111) surface parallel the decagonal surface 
of the quasicrystaline substrate, and whose azimuthal orientations in equal in- 
crements of 36° results in a pseudodecagonal symmetry (FCC domain phase) 
[S] [TD]. Our simulations show that the QC-like phase emerges at fast thermal 
relaxation times and weak adatom interactions, with the growth mode being 
layer by-layer. The FCC domain phase occurs when either the relaxation times 
or the relative strength of the adatom interactions is increased along with a 
transition from layer by layer growth to a cluster type growth that seems to 
depend predominantly on the thermal relaxation times. 

We are thus able to reproduce numerically the formation and size selection 
of domains, their distinct orientations and have shown how the depositional 
dynamics of adatoms on the substrate determines its morphology. Our ex- 
perimental and computational efforts in determining the structure of growing 
adlayers have shown that there is a structural registry between the crystal and 
the quasicrystal on a local scale. Experimentally observed diffraction patterns 
contain an orientational smearing of the structure with a well-defined length 
scale associated with adatom spacing and symmetry. 

The results indicate that epitaxial oxygen atoms chemisorbed on the qua- 
sicrystalline surface experience competing interactions. While they favor or- 
dering in a stable hexagonal mesh, the quasicrystalline surface template they 
condense on, force them into an aperiodic order. The system finds the best 
compromise by partially satisfying both conditions and breaks up into domains 
ordered in a hexagonal structure, where each domain remains locally commen- 
surate and in registry with the substrate. This registry results in five distinct 
orientations for the oxygen nanocrystals. The size of these domains is deter- 
mined by the interfacial strain energy, as this size increases the island edges 
get rapidly out of registry with the substrate and it becomes energetically more 
favorable to break-up into locally commensurate domains. 

In heteroepitaxial systems the interfacial strain energy increases with the 
film thickness, and beyond a critical thickness strain built in the film causes 
the film to relax to its stable bulk phase by creating misfit dislocations at the 
interface. In monoatomic systems on decagonal quasicrystalline surfaces, this 
critical thickness does not seem to extend beyond a monolayer thickness, and 
the system relaxes into large hexagonal domains. It is likely, however that. 
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in some carefully chosen binary and ternary alloy systems the critical thick- 
ness may increase significantly, and one may be able to grow epitaxy stabilized 
quasicrystalline phases, that do not exist in nature. 

Computations were mainly done using the Gilgamesh and Kassandra com- 
puter clusters at the Feza Giirsey Institute and Bogazigi University. This work 
has been funded in part by grant 08B302 of Bogazigi University. Financial 
support by Schweizerischer Nationalfonds is gratefully appreciated. 
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